Nonlinear smectic elasticity of helical state in cholesterics and helimagnets 
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General symmetry arguments, dating back to de Gennes dictate that at scales longer than the 
pitch, the low-energy elasticity of a chiral nematic hquid crystal (cholesteric) and of a Dzyaloshinskh- 
Morya (DM) spiral state in a helimagnet with negligible crystal symmetry fields (e.g., MnSi, FeGe) 
is identical to that of a smectic hquid crystal, thereby inheriting its rich phenomenology. Start- 
ing with a chiral Frank free-energy (exchange and DM interactions of a helimagnet) we present a 
transparent derivation of the fully nonlinear Goldstone mode elasticity, which involves an analog 
of the Anderson-Higgs mechanism that locks the spiral orthonormal (director/magnetic moment) 
frame to the cholesteric (helical) layers. This shows explicitly the reduction of three orientational 
modes of a cholesteric down to a single phonon Goldstone mode that emerges on scales longer than 
the pitch. At a harmonic level our result reduces to that derived many years ago by Lubensky and 
collaborators [J . 
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Dating back to the original cholesterol liquid crystal 
discovered by Reinitzer, chirality plays a central role in 
modern study of liquid crystals [2]. It is equally impor- 
tant for understanding helical states of noncentrosym- 
metric magnets (e.g., MnSi)[!, 0], driven by a chiral 
Dzyaloshinskii-Morya (DM) interaction [|| 

Among a wealth of induced phenomena^ chirality 
converts uniform nematic and ferromagnetic phases into 
states in which the orientational (nematic or spin) field 
twists periodically, thereby leading to a variety of spa- 
tially modulated phases, such as the cholesteric, two- 
and three-dimensional "Blue" phases 0,1(3,0] and recently 
discovered Skyrmion line crystals in MnSi helimagnet [8] . 
These spontaneously break the translational and rota- 
tional symmetries, forming liquid-crystalline structures 
that are periodic along 1, 2 or 3 dimensions. A cholesteric 
and helical phases are the most ubiquitous of these, char- 
acterized by a biaxial order with 



h (r) — x\ cos(q • r) + i 2 sin(q • r) 



(1) 



breaking the translational symmetry along a single, spon- 
taneously selected axis, £3 = z, where x\ , x% , xs = x± x x^ 
form an orthonormal triad that is constant in the ground 
state. 

General symmetry arguments Q and an explicit deriva- 
tion at a harmonic level [ij applied to this spontaneously 
layered helical state predict that at scales longer than 
the helical pitch, the low-energy (Goldstone mode) elas- 
ticity is identical to that of a smectic liquid crystal [2]. 
The three orientational degrees of freedom defining the 
cholesteric at short scales thereby reduce to a single 
smectic-like phonon mode, that emerges on scales longer 
than the pitch. The associated enhanced fluctuations 
in[l, [To| and disordering of these states leads to rich 
phenomenology 0, H2, an d m the case of MnSi is be- 
lieved to be associated with the striking observation of 
the non- Fermi liquid behavior [ll). 



Although by now quite familiar, the symmetry break- 
ing of the helical state falls outside the conventional G/H 
paradigm [l2j. Despite fully breaking the 0(3) rotational 
symmetry of the Euclidean group the state is character- 
ized by only a single U{1) Goldstone mode, x, t ne s pi~ 
ral's phase related to the smectic phonon u = — x/Qo- 
As we will show below, the absence of the two ad- 
ditional orientational Goldstone modes is best under- 
stood as a mathematical equivalent of the Anderson- 
Higgs mechanism [2I [T2I [l3j that gaps them out. 

In this Letter we explicitly show how this single low- 
energy helical mode emerges and derive its nonlinear 
smectic energy functional, expected from the underly- 
ing rotational symmetry 0, [Hj]. Because the latter leads 
to harmonic phonon fluctuations that diverge in three 
dimensions and below, the inclusion of nonlinearities is 
essential for a sensible and self-consistent description, as 
anticipated long ago by Lubensky, et al.,[l[ and Grin- 
stein and Pelcovits[9lll0l]. With the neglect of crystalline 
anisotropy all of our cholesteric results apply equally well 
to the description of the low-energy bosonic modes of the 
helical state in the DM helimagnets such as MnSi@,||. 

The helical texture no( r ) minimizes the chiral Frank- 
Oseen free-energy density of a chiral nematicjl, HH, [H| 



H = 



^K S (V -h) 2 + ^K b (hx Vxn) 2 , 
+ -K t (h- V x n + q ) 2 , 



(2a) 



= -K[(d in3 ) 2 +2q h-V xn + q 2 ] , (2b) 

where in the second form we focused on the isotropic 
limit, K s = Kb = K t = K and dropped the total deriva- 
tive saddle-splay (Gaussian curvature) contribution, that 
reduce the Hamiltonian to that of DM ferromagnet in the 
absence of crystal symmetry breaking fields Within 
this approximation the space-spin (r — n) coupling only 
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enters through the chiral (second) term, with the elastic 
(first) piece explicitly exhibiting independent rotational 
invariances of space, r and of the director, n (spin in the 
MnSi magnet context). 

We now look at long- wavelength, low-energy Goldstone 
modes excitations about the helical ground state, no(r). 
A general state is described by 

ri(r) = ei(r) cos(q ■ r + x(r)) + e 2 (r) sin(q • r + x(r)), 

(3) 

where ei(r), e 2 (r), es(r) = ei(r) x e 2 (r) now constitute 
a spatially dependent orthonormal frame, describing the 
orientation of the local director (spin) helical plane, that 
is independent of the helical axis set by q • r + x( r ) — 
const. The phase x(r) = — qou(r) of the chiral helix 
also defines the phonon field u(r) of the helical layers, 
which on scales longer than the pitch a = 27r/go de- 
fine the smectic-like displacement of these helical phase 
fronts. Thus, altogether on the intermediate scales there 
are three independent orientational degrees of freedom, 
x(r) and ea(r). The azimuthal angle <fi(r), defining the 
orientation of the e% t 2 around e 3 is redundant to x( r ); as 
it can be eliminated in favor of it via a local gauge-like 



transformation on ei f2 . Although naively, the low-energy 
coset space is isomorphic to (Si x S 2 )/Z 2 — RP 3l as we 
will see below, only a single Goldstone mode, x( r ) wm 
survive this helical symmetry breaking. 

Substituting n(r) from Eq. © into the free-energy of 
the chiral nematic (helimagnet), Eq. (I2bp . and using 

dihj = (qi +9ix)[-e y sin(q-r + x) + e 2 j cos(q • r + x)] 
+d l ei J cos(q • r + x) + d l e 2j sin(q ■ r + x), (4) 

together with 



d l eij = aiB 2 j + cuesj , (5a) 

9ie 2 j = -a t e 2j + c 2i e 3j , (5b) 

and gauge field "spin-connections" 

a>i = e 2 • diSi, (6a) 

Cu = -ei ■ <%e 3 , (6b) 

c 2 i = -e 2 • dii 3 , (6c) 



we find 



H = 



K 

y 



(aie 2j + Ciezj) cos(q • r + x) + (-adij + c 2i e 3 j) sin(q • r + x) 



+ diX) [ - eij sin(q • r + x) + e 2j cos(q ■ r + x)] 

+go(e 2l e 3j - e 3l e 2j ) cos(q • r + x) + qo(hi£ij ~ eue 3j ) sin(q • r + x) 

K K K 

— (Vx + a + q - q e 3 ) 2 + — (ci + q i 2 ) 2 + — (c 2 - q ei) 2 



(7a) 
(7b) 
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In obtaining Eq. (|7b[) , we dropped the constant and oscil- 
latory parts, that average away upon spatial integration. 

We first note that the requirement of well-defined he- 
lical phase fronts, i.e., the absence of dislocations and 
disclinations in the layer structure, can be enforced by 
the compatibility condition V x a = 0, consistent with 
Mermin-Ho relation [15J. This allows us to take a = V</> 
and thereby shift ("gauge") away </>(r) in favor of x( r )> 
according to x( r ) + 4>( r ) ~^ x( r )- 

Without loss of generality, we next take q = qoz with 
z defining the orientation of the helical axis in the lab- 
oratory coordinate system x,y,z. The long-wavelength 
free-energy density, reexpressed in terms of the smectic- 
like phonon field u(r) and the local nematic helical frame 
orientation e 3 then reduces to 



2 

, Kq 



(Wu + e 3 - z) 2 + ^(cl + c 2 2 ) 



(ci • e 2 - c 2 • ei). 



(8) 



Clearly the first term, above, accounts for the energetic 
cost of the deviation of the local nematic frame e 3 from 
the local orientation of the helical layers. A minimiza- 
tion of this term (or equivalently at long wavelengths, in 
a statistical mechanical treatment integrating out the in- 
dependent e 3 (r) degree of freedom) , at low-energies locks 
the orientations of the cholesteric layers and the nematic 
frame. In a perturbative treatment this leads to the ex- 
pected relation 



u « -e 3J _, 



(9) 



that is an example of a Higg's-like mechanism (akin to 
thermotropic smectic liquid crystals 0, [H, HH), that at 
long scales effectively gaps out the orientational Gold- 
stone modes. 

The exact minimization over the unit vector e 3 (r) can 
also be carried out using a Lagrange multiplier A to im- 
pose the constraint £3-63 = 1. Minimization over e 3 (r) 
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gives 

Vw + e 3 - z = -Ae 3 , (10) 

with the solution 

A + l = ^{Vu - zf = VI - 2u zz , (11a) 

e 3 = (z-Vu)/y/l-2u zz , (lib) 

where 

u zz = d z u~\{Vu) 2 (12) 



is the standard fully nonlinear smectic strain tensor, 
that encodes the full rotational invariance of the heli- 
cal stated [12]. Using Eqs.CU} to eliminate A and es(r) 
in favor u zz (valid at long scales), we find 



n = 



Kql 



K 



(VI - 2u zz - l) 2 + — [(<§! • die 3 )(ei ■ 9;e 3 ) + (e 2 • die 3 )(e 2 ■ ^e 3 )] + —^-eue 2 j(didjU - djdiu), (13a) 



ju 2 zz + ^-[didfuf, 

|[a zM -i(v u ) 2 ] 2 + |(vi u ) 2 



(13b) 
(13c) 



where to obtain our main result, Eq. (|13c[) we used the 
condition of well-defined helical layers with no disloca- 
tions in u(r), i.e., a single-valued phase field x( r )> ne ~ 
glcctcd the boundary terms, expanded to lowest order in 
the nonlinear strain tensor u zz , and defined the compres- 
sional and bending elastic moduli 



B 



Kql K = K/2. 



(14) 



Thus, as advertised, we have demonstrated that on 
scales longer than the helical pitch 2ir/qo, the low-energy 
deformations (Goldstone modes) of the helical state are 
characterized by a fully rotationally invariant, nonlinear 
smectic elastic theory (3, ||| [l2j]. The latter can be de- 
rived by spontaneously ordering the density p(r) of the 
isotropic fluid into a one-dimensional periodically mod- 
ulated state (smectic), characterized by p(r) = p + 
p q cos[(q • r — <7ou)][2(- Alternatively, the above nonlinear 
compressional form (first term in (|13a[l ) emerges directly 
from the de Gennes' gauge theory of the smectic[2[, after 
condensing p q = jpqje 1 *^^ to give 



n sm = -B(V$ - nf 



(15) 



and then minimizing over the unit director n(r) as in 

(nni-dni). 

We note that as required, T-L in Eq. (113a[) is a function 
of the fully nonlinear strain u zz as it must to preserve 
full rotational invariance. Furthermore, it is a nonlinear 



function of this strain, that reduces to the familiar "har- 
monic nonlinear" formQ in Eq. (|13c[) only for small u zz . 
It is worth observing that through the introduction of 
the phase field $ = — z + w, this nonlinear in u zz term in 
(|13ap can be written as (IV^I — l) 2 . This form has been 
used in recent analyses of various nonlinear properties of 
smectic. fl7jj 

Our result in Eq. (|13cp then in turn implies that 
the helical state inherits all the novel nonlinear elas- 
tic effects previously discovered in the context of con- 
ventional, thermotropics and lyotropic smectic liquid 
crystals and other spontaneously layered states that 
emerge from an isotropic state. These include thermal 
fluctuations d, [l(| [3 and heterogeneous fl9i anomalous 
elasticity effects, the undulation instability [20], and many 
others. 

One important distinction from a conventional smec- 
tic, however, is the underlying chirality of the helical lay- 
ered state. Although as in the chiral smectic the effective 
Anderson-Higgs mechanism expels the expression of chi- 
rality (e.g., twist) inside the helical state (as the magnetic 
flux density [twist of the vector potential] is expelled from 
the Meissner state) ddHHij], inclusion of the chiral terms 
(encoded in the departure from the Vxa = condition 
used to get to ([5])) is essential to understandin g th e topo- 
logical defects and melting of the helical state 21 [ . 
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